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Abstract 

The class HCM consists of all non-negative functions / on (0, oo) 
such that f{uv) ■ f{u/v) is completely monotone with respect to (wrt) 
w = V -\- 1/u, for every fixed m > 0, and has been extensively studied 
for a long time. It is closed wrt a number of useful operations, such 
as pointwise limits and products. It is also closed wrt products and 
quotients of independent random variables, some changes of variables 
and the Laplace transform. We consider its multivariate (bivariate) 
counterparts MVHCM (BVHCM) and study some of their properties. 

In particular, we show that MVHCM is closed wrt the Laplace trans¬ 
form and use this to define the class BVGGC-L of bivariate random 
vectors with Laplace transform in BVHCM. Then BVGGC-L contains 
Bondesson’s class of random vectors in BVGGC in the strong sense 
and is contained in the corresponding class BVGGC in the weak sense. 

We also show that BVHCM, in contrast to HCM, is not closed wrt 
multiplication of independent bivariate random vectors. 
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1 Introduction. 

We consider various classes of multivariate probability densities f{x) on the 
positive cone = (0, oo)” in R”. In particular, we consider the class of 
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completely monotone functions (CM) and the class of hyperbolically com¬ 
pletely monotone functions (HCM), introduced in [5], Ch. 5. A function / 
is completely monotone (CM) if / is on and (—l)l“l • D°‘f{x) > 0, 
for all multiindices a = (ai,... ,0^) of nonnegative integers ai, 1 < i < n, 
and |q;| = ai a„. The CM-functions on (0, oo) were characterized 

by Bernstein [2]) P- 56 as Laplace transforms of nonnegative measures. See 
also Widder m Theorem 12 a,b and Bochner [3] Theorem 4.2.1. The class 
of CM-functions and the related classes of Stieltjes and Bernstein functions 
are studied in [9] and m- 

The class of hyperbolically completely monotone (HCM) functions on (0, oo) 
was introduced in Bondesson [5] for the study of generalized gamma con¬ 
volutions (GGC). A function f(x) defined on (0, oo) belongs to HCM if 
f{uv)-f{u/v) is CM inw = v+l/v, for all fixed n > 0. We write X ~ HCM 
for a random variable A, if its density function f{x) is HCM. The class HCM 
is extensively studied and characterized in [5] (and there denoted by C). 

The bivariate analogue BVHCM of HCM was introduced in [1], p. 193 (and 
there denoted by C 2 ). In contrast to the univariate case, the class BVHCM 
is much less studied. A comprehensive study of continuous bivariate distri¬ 
butions, without mentioning hyperbolicity, is found in [T]. 

The class MVHCM of n-variate hyperbolically completely monotone func¬ 
tions f{x) = f{xi,... ,Xn) defined in i?” and the corresponding class of 
random vectors is defined in analogy with the bivariate case (see Section 
2). We state and prove some of its properties related to marginal and con¬ 
ditional distributions. It was conjectured by Bondesson [6], p. 3764, that 
MVHCM is closed wrt the Laplace transform. We conhrm the conjecture 
in Theorem 1 and use this fact to define the class BVGGC-L of bivariate 
densitiy functions with Laplace transform in BVHCM. Then BVHCM-L 
contains Bondesson’s class of random vectors X ~ BVGGC in the strong 
sense (Theorem 2) and is contained in the weak class BVHCM. 

It is well known that if A ^ HCM and Y ^ HCM are independent ran¬ 
dom variables, then XY ~ HCM and X/Y ~ HCM, [5] Theorem 5.1.1. 
The analogous result for independent random vectors (A, V) ~ BVHCM 
and {X',Y') ~ BVHCM was posed as an open question in [6], p. 3764. 
Our main result answers this question in the negative (Theorem 3). The 
reason is that if f{x,y) denotes the density function of the product vector 
{XX',YY'), then f{vi,V 2 ) ■ f{l/vi,l/v 2 ) depends in a crusial way on the 
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quadratic term v\lv 2 + V 2 lv\ = wi-w^ — W 2 , where as usual wi = vi + l/vi, 
W 2 = V 2 + I/v 2 and tea = ui /u 2 + ^ 2/^1 . We give an explicit example that 
exploits this fact. 

2 Preliminaries. 

Let BP denote the n—dimensional Euclidean space and consider non¬ 
negative functions f{x) = /(xi,... ,Xn) defined in the positive cone BP\_ = 
{x G R'^;xi > 0, 1 < f < n} of R^. Integration wrt Lebesgue measure in R^ 
is written f /(x) dx. Our notation for measures and integrations are stan¬ 
dard and for the necessary background in probability theory we rely on [5] 
and [9]. The following definition of the class of the n-variate hyperbolically 
completely monotone functions is in analogy with the bivariate case in [7]. 

Definition 1 A n-variate density f defined on Rfi is MVHCM (MVnHCM) 

if 

fiuiVi,U2V2, . . ■,UnVn) ' f {ui f Vi, U 2 IV2, ■ ■ .,UnlVn) (1) 

is completely monotone as a function of Wi = Vi + 1/vi, 1 < i < n and 
Wij = Vijvj + Vj jvi, I < i < j < n, for all positive numbers Ui, 1 < i < n. 
In the bivariate case n = 2 we denote wi ^2 by w^. 

As an example of functions in MVHCM we consider 

n 

/(3;) = ■ (1 ai,iXl H-h an,iXn )~'^', 

1 

where Oi > 0, a* > 0, 1 < f < n and 7 > 0, properly normalized. Then the 
product in ( 1 ) takes the form 

n n 

^ . y,. ^ aiOjUiUj ■ (2) 

2 = 1 2=1 i<j 

which clearly is MVHCM, and can be represented as 

/ n 

exp ( - Ai • rcj - y] XijWij) dn{Xi, A^), 

1 i<j 

for some nonnegative measure n depending on ai,Ui,ai,'y and m = ( 2 ). 
This can be seen directly or by [3] Theorem 4.2.1. 
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3 Main results. 


We begin this section with a number of properties of random vectors {X, Y) 
with densities / in BVHCM. The first result states that the conditional 
density fx\Y=y{x), the marginal density fx{x), and the quotient density 
fx/vi^) 3-11 HCM and that BVHVM is closed wrt taking q-th powers, 
for \q\ > 1. 

Proposition 3.1 (a) Let {X^Y) ~ BVHCM and denote the density func¬ 
tion by f. Then fx\Y=y{x), fx{x) and fx/xi^) ^xe HCM. 

(b) Let {Xi,X 2 ) ~ BVHCM, then {Xf^,X:^^) ~ BVHCM. 

(c) Let {Xi,X 2 ) ~ BVHCM, then (Xf,X|) ~ BVHCM for |g| > 1. 

Proof. To prove (a) we first let Fi{x) = fx\Y=yix) = f{x,y), for some hxed 
y > 0. Then Fi{uv) ■ Fi{u/v) = f{uv, y ■ 1) ■ f{u/v, y/1) is CM in n + l/v, 

OO 

since / is BVHCM. Next let F 2 {x) = fx{x) = f f{x,y) dy, then 

0 


F 2 {uv)F 2 {^) = J J f{uv,y)f{^,z) dydz = 2j J f{uv,st)f{^,^)-^dsdt, 
0 0 0 0 

by Fubini’s Theorem and the hyperbolic change of variables y = st, z = s/t 
with Jacobian —2sft. For hxed u and s there exists a nonnegative measure 
n, depending on u and s, such that the last integral equals 


2 - J sds j y j exp ( - Ai(n + i) - A2(t + ^) - A3(^ + ;^)) di/(Ai, A2, A3), 
0 0 

since / is BVHCM. The last two terms in the integrand of the inner most 
integral, with reversed signs, becomes 

t • (^2 H-) + — • (A2 + A3n) = p ■ (A2 + A| + A2A3 • {v H—)) 

V t V 

after a change of variables t = p-{X 2 + X 3 v). Putting the pieces together, this 
proves that F 2 is CM in the variable v + l/v. Finally, let F^{z) = fx/xi^) = 
/o°° y ■ fizy, y) dy, then F 3 {uv)F 3 {u/v) equals 


OO OO 



0 0 


dxdy xy ■ f{uvy, y) ■ f{-x, x) 


f f , dt , „,u s s, 

= 2-y j s^ ds--f{uvst,st)-f{---,j), 

0 0 
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by Fubini’s Theorem and the same hyperbolic change of variables. The last 
two factors are CM in the variables vt + l/vt, t + l/t and v + l/v and the 
proof can be completed in analogy with the case F 2 . 

For (b), we let / denote the density function of [Xi,X 2 ) and put Z = 
Then fz{xi,X 2 ) = and / HCM implies 

that fz is HCM. 


In view of (b), it is enough to prove (c) for q > 1. Let / be as in (b) and put 
Z = X|), then/^(xi, X 2 ) = ^x^ 2 '^ ^ •/(x}^'^, x^^'^). The expression 


fziuiVi,U2V2) ■ /(—, —) 

Vi V2 




fz{u{‘v^ , ^ 
can, for fixed ui and U 2 -, be represented by 


■ fz{ui 




Uo 


0 


j exp [ - Ai • {v^ ") - A 2 • (u| +X 2 - A 3 • ((^)'' + (—)")] 


for some nonnegative measure dz^(Ai, A 2 , A 3 ). The terms is the exponent, 
with reversed sign, are Bernstein function in the variables ui + — ,U 2 + — 
and ^ respectively, by [3j, p. 60. A 

Remark 1. The statements in (a) about fx\Y=yix) and fx{x) also hold when 
the random variables X and Y are replaced by random vectors. The details 
are left to the reader. We show by an example that fx\Y=y{x)-, fY\x=xix), 
fx{x) and fY{x) HCM does not imply that f(x,Y) BVHCM. 

Example. Define f{x,y) = c ■ {1 + x + y + kxy)~^ , where 7 > 0, c is a 
normalizing constant and A; > 1. Clearly, fx\Y=y{x) and fxix) are HCM. 
Put H = f{vi,V 2 ) ■ f{l/vi, I/U 2 ), then 

FI ~ (3 + + (1 + k){wi + W2) + W3 + k{wi ■ W2 - W3)) ^ 

and > 0. Thus / is not BVHCM, which completes the example. A 


The following multiplication theorem seems to be new and answers a ques¬ 
tion posed by Bondesson (personal communication). 

Proposition 3.2 Let (Vi,X 2 ) ~ BVHCM and Y ~ HCM, then Z = 
{Y -XuY ■X 2 ) BVHCM. 

Proof. Denote the density functions of (Vi, V 2 ) and Y by f{x,y) and g{x), 
respectively. Then 

00 

fzix,y) = J •5(A)^ 

0 
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and fz{uiVi,U 2 V 2 ) ■ fz{uilvi,U 2 lv 2 ) becomes 


J = 2 


OO OO 


0 0 


UlVl U2V2-. 
xy ' xy ' 


/( 


uiy U2y ^ 

ViX ’ V 2 X ‘ 


X dx dy 
9 {xy)-g[-) -3 — 
y x-^ y 


( 3 ) 


after a hyperbolic change of variables. By our assumptions on / and g, 
this integrand is, for fixed mi,U 2 and x, CM in the variables vi/y + y/vi, 
V2/y + y/v2-, vi/v2 +V2/V1 and y+ l/y axid can be represented as f e~^ dv, 
where 


y Vi y V2 V 2 Vi y 1 

and is a nonnegative measure. After a change of the order of integration 

OO 

in (3), it is enough to prove that the inner most integral f e' -E dy ig CM. 

0 ^ 

Three of the terms in E contain the variable y and can be written 

y ■ (-1-h A4) H— • (AiUi + A2U2 + A4) = 

vi V 2 y 


— P ' (^1 T A 2 + A 4 + A 1 A 2 • 1^3 + A 1 A 4 • W\ + A 2 A 4 • W 2 ) + 1/p, 

after the change of variables y = p - (AiUi + A 2 U 2 + A 4 ), from which it is clear 
that J is CM and fz is BVHCM. A 


Our next result states that the class MVHCM is closed wrt the Laplace 
transform and answers a question in [^, p. 3764. In the univariate case this 
follows from the multiplication theorem for independent random variables 
with densities in HCM, see [5] Theorem 5.1.1. 

Theorem 1 Assume that f is MVHCM, then the Laplace transform of f 
is also MVHCM. 

Proof. We give the proof for n = 2 only, the general case being similar. Let 
/ be as in the theorem, denote the Laplace transform of / by T and define 
i7 = F(uiUi,U 2 U 2 )-F(|^,^). Then 

H= dsidS2dtidt2 eX.p{-SiUiVi-S2U2V2-ti — -t2 — )-f{si,S2)-f{tl,t2), 
J Vi V2 

where the integration is over R'^. Now we make the hyperbolic changes of 
variables si = xy,ti = x/y and S 2 = zw,t 2 = zjw with Jacobians —2x/y 
and —Izjw respectively and get 

H= [■■■ •/(-,-), 

.J J y w 


6 




which is our starting point. For fixed ui,U 2 -,x and z we get by Bernstein’s 
Theorem, [3] Theorem 4.2.1, that the last two factors can be represented as 


f{xy,zw) ■ 

y w 


J Q->^v{y+i-/y)-^2-{w+l/w)-X3-{y/w+w/y) A 2 A 3 ) 


for some nonnegative measure v depending on tti, U 2 ) 3 : and z. Changing the 
order of integration we find that it is sufficient to prove that the inner most 
integral in H 

j_ f dydw ^-uix{viy+:;^)-U2z{v2W+:;^)-\i-{y+^^)-\2-{w+^)-X3-{^^ + ^) 

J yw 


is completely monotone as a function of rci = ui + 1 /ui, W 2 = U 2 + l/u 2 
and W 3 = vi/v 2 + V 2 /vi. The exponent in the integrand of I can be written 
(with the sign changed and new names on the variables) 


J = A{viy + ^—) + B{v 2 W H—3—) + C{y+-) + D{w + —) + E{— + —) 
viy V 2 W y w w y 


We are going to perform two changes of variables in J and start with the 
y—variabel. The six terms in J that contain a factor y or 1/y can be written 

Ji = y{Avi + C + —) + -(— + C + Ew) 
w y vi 


and we make the change of variables y = p ■ {A/vi + C + Ew). Then 
dt/t = dp/p and 


j^ = p.(A^ + C‘^ + E‘^ + AC ■wi + CE-{w + -) + AE- (wvi + —)) + -. 

W WVl p 

We treat the remaining four terms in K together with the terms in Ji 
containing a factor w or 1 /w 


J2 = w ■ {Bv2 + {D + pCE) + pAEvi) H-• (- {D + pCE) + pAE —) 

W V 2 Vl 

in a similar way. The change of variables w = 5 ■ + {D + pCE) + 

pAE^), d 6/6 = dw/w, expresses J2 as a linear combination of wi, W2 and 
w^. Putting the expressions for Ji and J2 together we find that 


J = a + b ■ wi + c ■ W 2 + d ■ W 3 , 

where a,b,c,d are polynomials in Ai, A 2 , A 3 , ui.n 2 , x, z, p, 1/p, 5 and 1/5 with 
positive coefficients. It follows that /, and thereby also H, is completely 
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monotone in the variables r(;i, W 2 and w^. This completes the proof. 


A 


There is a close connection between HCM and the class of generalized gamma 
convolutions (GGC) introduced by Thorin [11], [T2| in 1977. The class 
GGC consists of the limit distributions of sums of independent Gamma 
distributed random variables, see [5|. Bondesson proved that a density f is 
GGC if and only its Laplace transform is HCM, |3] Theorem 5.3.1. Two 
classes of bivariate GGC were defined in [5], pp. 46 - 47. A distribution 
on R'^ is BVGCC in the strong sense if it is the limit distribution for a 
sequence of random vectors ^1=1 '''^here Yj are 

independent, unit scale Gamma variables and cij-, C 2 j > 0. The distribution 
of a random vector (Ai, A 2 ), where Ai, A 2 > 0, is BVGGC in the weak sense 
if Cl ■ Ai + C 2 • A 2 ~ GGC for all ci, C 2 > 0. In view of Theorem 2 and the 
characterization of GGC mentioned above, we define BVGGC-L as the class 
of bivariate densities on for which the Laplace transform is BVHCM. 
Then BVGGC-L is closed wrt addition of random variables and we have the 
following result. 

Theorem 2 If X ~ BVGGC in the strong sense, then also X ~ BVGGC- 
L. 

Proof. Let be the bivariate random vector in the definition of BVGGC. 
Then the Laplace transform of at the point (si, S 2 ) is 



which is BVHCM, since BVHCM is closed wrt multiplication of functions. 


Taking limits as n —)• 00 completes the proof. 


A 


It is easy to see that A ~ BVGGC — L implies that A ~ BVGGC in the 
weak sense, which gives the inclusions 


BVGGC strong C BVGGC - L c BVGGC^eak. 


The multiplication theorem for HCM states that if A ~ LICM and Y ~ 
HCM are independent random variables, then also A • V ~ HCM and 
X/Y ~ HCM, [ 5 ] Theorem 5.1.1. The corresponding result for bivariate 
random vectors was posed as an open question in [6] p. 3764, but so far no 
proof has appeared. We show by an example that the such a statement is 
in general false. 
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Theorem 3 There exist independent random vectors {X,Y) ~ BVHCM 
and {X',Y') ~ BVHCM such that the product {X ■ X' ,Y -Y') does not have 
a density function in BVHCM. In particular, this is the case if {X,Y) and 
{X',Y') have densities f{x,y) = g{x,y) = y~^ ■ ^ respectively, for 

k > 0 sufficiently large. 


Proof. We set the stage by letting / and g denote the density functions for 
the two arbitrary independent random vectors {X,Y) and {X',Y'), respec¬ 
tively, on Since X and X' are independent we get 


OO OO OO 

fxx'{x)= [ — [ f{x/s,f) [ g{s,u)dtdu 


0 


0 


0 


and similarity for fvY'ix), By the independence assumption, we find the 
density function for the random vector {XX' ,YY') to be 


F{x, y)= j ^ j y/(x/s, y/t) ■ g{s, t). (4) 

0 0 


Now we put f{x, y) = g{x, y) = y ^ • e * and get 


OO OO OO OO 



0 0 0 0 


From the last formula it can be shown that F{x, y) is finite for all x,y > 0. 
Let vi,V 2 >0 and define J = F{vi,V 2 ) ■ F{l/vi, I/U 2 ), then 


J = 



OO 


/ 


dt - 
—e 
t 


•"1 j. '" 1 * 
s SV2 


-s-kf 



dv — 1 — 

_g V\U V\U V 

V 


We will show that J is not CM in the variables rci = vi + l/vi,W 2 = V 2 + I/V 2 
and tC 3 = V 1 /V 2 +V 2 /V 1 , provided k is large enough. As a first step we express 
J as a function of wi , W 2 and wj, by two hyperbolic changes of variables 
followed by two changes of variables similar to the ones in the proofs of 
Proposition 3.1 and Theorem 1. We begin by putting s = xy,u = x/y and 
t = zw,v = z/w, with Jacobians —‘Ixjy and —^zjw respectively, which 
after some algebra gives J = 4 ■ f f f f e~^dxdydzdw/xyzw and 


x y vi y x V 2 y 


V 2 V, kx ,v w, 

-) + — •(- + -) 
viw z w y 


9 



,1 , Z V 2 1/1 , z Viw , X , 

= y • (—"Wi + a: + k - h K - j H — • (— + x + k - h k—w). 

X X viw z w y X X V 2 z 

We make a change of variables by putting the second term equal to l/p and 
get 



+ kz{^^ + + k—{w + —) + rci + kw-i) + l/p. (5) 

V\W V 2 z w 

The four terms in (4) that contain the variable w are now treated in the 
same way. These terms are p times 


,kz 1 k 1 , ui kx‘^ ^ 1 ,kz k , uo kx'^^ 

+-h kz— +-) + - • (^V2 + -vi + kz— +-). 

X^ V2 Z Vi V2 Z ' W ^X^ Z Ul Z ' 


A change of variables by putting the last term equal to 1/(5 gives 


E' = (5( 


k'^z^ k"^ 


+ 


k‘^x 


H—^ k z + (— K — h 




) • rci + (A:^x^ + • W3+ 


+k‘^ ■ W 2 + k‘^ ■ {— + ^)) + 1/(5. 
V 2 v( 


The crusial term here is the nonlinear term ^ + p- = wiws — W 2 , which 
turns out to be responsible for J not being CM. It will be our main concern 
in the rest of the proof. We also note that W 2 is cancelled and that the terms 
that contain wi and ws are of two types. The first type is terms with no x— 
and z—variables, they are p ■ wi + pk ■ W 3 + pdk"^ ■ W 1 W 3 . The second type of 
terms are linear combinations of wi and W 3 with coefficients depending on 
p, 5, X and z. If we put all this together we get the following expression for 
E 


E = gi{x,z,p,5) + g 2 {wi,W 3 ,x,z,p,S) + p ■ wi + pk ■ W 3 + pdk"^ ■ W 1 W 3 , 

where gi and g 2 are polynomials in the indicated variables, 1/p and 1/5. 
Carrying out the x— and z—integration gives the formula 

00 00 

J = j J — y rua, p, 5) ■ ( 5 ) 

0 O 

on which the remaining part of the proof is based. A careful study of the 
expression for E shows that J and all its derivatives are finite for wi,W 3 > 0. 
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We claim that the second order mixed derivative of J is negative for wi , W 2 
close to the origin and k sufficiently large, which proves that J is not CM 
and F is not BVHCM. 

In the following we denote the integration in (5) by f dm and suppress the 
variables p, 5 in G. Denoting the derivative of J and G wrt wi by Ji and 
Gi, respectively, gives 

Ji = J . (^Gi{wi,W3) - G{wi,W3) ■ {p + p6k‘^W3)) 

and, by letting wi 0, 

Ji = J ■ {Gi{0,W3) - G{0,W3) ■ {p + pSk^ws)). 

We let Ji ^3 denote the derivative of Ji wrt W 3 , and similarity for G, and 
suppress the tci-variable from the notation. Then 

•^1,3 = j ■ [Gp3(u)3) - G{w 3 ) ■ pdk"^ - Gsiws) ■ {p + pdk'^ws) 

-pk ■ {Gi{w3) - G{w3) ■ {p + pSk'^W 3 ))]. 

Now we let /c —7> 00 and tcs ^ 0 in such a way that k^ws —>■ 0. Then Ji ^3 is 
asymptotically 

J dm (Gi, 3 ( 0 ) - G(0) • p5k^ - G 3 ( 0 ) - pk ■ Gi(0) + G(0) • p^k) < 0, 

if k is sufficiently large, since all integrals are finite. This proves our claim 
and completes the proof of Theorem 3. A 


Remark 2. Let {X,Y) and {X',Y') be independent random vectors in 
BVHCM with density functions f{x,y) and g{x,y), respectively, and de¬ 
fine F{x, y) by (4). Then, in a way similar to the proof of Theorem 3, it can 
be proved that J = F{uiVi,U 2 V 2 ) ■ F{ui/vi^U 2 /v 2 ) equals four times 


00 00 00 00 


/ dx k dy k dz k dw uivi 

X J y J Z J w X y’ 
0000 


U2W 

Z W ^ XVl ’ ZV2 


g{xy,zw)g{-,—] 
y w 


(7) 

By assumption, the first two factors and the last two factors in the inte¬ 
grand in ( 6 ) can be represented as Laplace transforms of nonnegative mea¬ 
sures fiiXi, A 2 , A 3 ) and r' 2 (A 4 , A 5 , Ae) respectively. Again, as in the proof of 
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Theorem 3, we can rewrite J as 


J = 4 


OO OO 

0 0 


(8) 


where a, b, c, d, e are polynomials in p, 5 and Ai,..., Ae- The crucial quantity 
here is the coefficient e, which is a product of p, S and four of the A* : s. If 
i^i and 1^2 have no mass on the set where e > 0, then J is CM. The converse 
is probably true, but seems hard to prove in the general case. It follows 
from (7) that, although J is not in general CM, J is CM in each variable 
separately. 


Acknowledgement. The author thanks Professor Lennart Bondesson for 
pointing out the problems treated in this paper and illuminating discussions. 
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